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ABSTRACT 

Finite- amplitude  standing  waves  contained  within  an  air 
filled  rigid-walled  tube  at  ambient  temperature  and  pressure 
were  investigated  experimentally.  The  pressure  waveform  in 
the  tube  was  analyzed  for  subharmonic  content  for  comparison 
with  several  existing  theoretical  models.  Due  to  limitations 
of  the  experimental  apparatus,  the  only  model  for  which  the 
predicted  threshold  for  subharmonic  generation  could  be 
exceeded  was  that  of  Coppens  which  predicts  subharmonic  gen¬ 
eration  when  the  strength  parameter  exceeds  2.0.  Strength 
parameter  is  defined  as  where  M2  is  the  Mach  number  of 

the  driving  frequency,  3  =  1.2  in  air,  and  is  the  quality 
factor  associated  with  the  subharmonic  resonance.  Strength 
parameters  up  to  2.89  were  investigated  and  no  threshold 
effect  was  observed. 
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I.  INTRODUCTION 


The  problem  of  large  amplitude  acoustic  waves  in  dissi¬ 
pative  media  has  been  extensively  studied.  It  is  well  known 
that  an  initially  sinusoidal  signal  of  large  amplitude  can 
distort  so  that  the  resulting  waveform  contains  frequency 
components  other  than  the  original  frequency.  Usually,  these 
additional  frequency  components  are  integer  multiples  (har¬ 
monics)  of  the  original  frequency.  However,  under  certain 
conditions,  components  at  frequencies  lower  than  the  original 
signal  frequency  may  be  generated.  These  lower  frequency 
components  are  known  as  subharmonics . 

The  objective  of  this  research  is  to  investigate  possible 
subharmonic  generation  with  finite- amplitude  standing  waves 
in  air  at  ambient  temperatures  in  a  rigid  walled  tube  for 
comparison  with  several  existing  theoretical  models. 
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II.  BACKGROUND  AND  THEORY 


In  a  nondissipative  medium  the  hydrodynamic  equations 
predict  that  an  elastic  wave  will  change  waveform  as  it 
propagates.  In  real  absorptive  media  experience  tells  us 
that  only  those  waves  of  relatively  high  amplitude  exper¬ 
ience  a  change  in  waveform.  Finite-amplitude  waves  are 
those  waves  which  undergo  harmonic  distortion  due  to  self¬ 
interaction  . 

The  study  of  intense  waves  began  in  1868  when  Kirchoff 
[1]  used  the  nonviscous  hydrodynamic  equations  to  predict 
a  change  in  waveform  for  a  traveling  plane  wave.  The  in¬ 
vestigations  of  finite-amplitude  effects  in  traveling  waves 
continued  [2,  3,  4,  5]  but  it  was  not  until  1954  that  Keller 
[6]  developed  solutions  for  standing  waves  excited  by  a  piston 
in  a  closed  tube.  Keller  assumed  a  nondissipative  medium 
which  resulted  in  infinite  amplitude  at  resonance;  however, 
his  results  were  still  useful  for  frequencies  close  to 
resonance.  In  1960,  Keck  and  Beyer  [7]  developed  a  non¬ 
linear  solution  for  progressive  waves  in  a  viscous  medium 
involving  terms  up  to  the  sixth  order  by  using  pertubation 
methods  and  Fourier  representation  of  the  waveform.  Coppens 
and  Sanders  [8]  extended  the  Keck-Beyer  pertubation  approach 
to  include  both  the  bulk  losses  and  wall  losses  as  predicted 
by  Rayleigh- Kirchoff  [9]  and  showed  excellent  agreement  with 
experiments  conducted  with  standing  waves  in  a  rigid-walled 
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tube  at  low  levels  of  nonlinear  interaction.  The  experiment¬ 
al  work  by  Beech  [10]  with  supporting  calculations  by  Ruff  [11] 
showed  that  at  high  excitation  levels  a  significant  difference 
developed  between  theory  and  experiment.  Winn  [12]  experi¬ 
mentally  demonstrated  that  the  Rayleigh-Kirchof f  loss  mech¬ 
anisms  failed  to  describe  the  phase  relationships  observed 
in  the  distorted  waveform.  Coppens  and  Sanders  [13]  revised 
their  model  to  utilize  empirically  determined  losses  and 
resonant  frequencies,  thus  extending  the  excellent  agreement 
between  theory  and  experiment  to  higher  excitation  levels. 

A  further  extension  of  their  model  to  two-  and  three-dimensions 
was  experimentally  investigated  by  Lane  [14]  and  found  to 
successfully  predict  the  major  features  of  the  harmonic  con¬ 
tent  for  finite  amplitude  standing  waves  in  a  rigid  walled 
cavity.  De  Vail  [15]  found  that  if  degeneracies  existed, 
the  model  failed  to  account  for  the  experimentally  observed 
excitation  of  non-family  modes.  The  model  was  again  revised 
and  experimentally  tested  by  Slocum  [16]  and  was  found  to 
accurately  predict  the  major  features  of  the  harmonic  content 
of  the  finite  amplitude  standing  wave  whether  or  not  degen¬ 
eracies  were  present. 

Investigation  into  the  nonlinear  generation  of  subharmonic 
frequencies  began  in  1965  when  Korpel  and  Adler  [17]  and 
Breazeale  and  McCluney  [18]  independently  observed  the  effects 
of  acoustic  subharmonic  generation  on  optical  diffraction 
patterns.  Cook  [19]  attempted  a  mathematical  explanation  of 
this  phenomenon  based  on  the  finite- amplitude  plane  wave 
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interaction  in  a  nondiss ipat ive  medium.  A  spectrum  analyzer 
was  used  to  investigate  subharmonic  generation  in  water  and 
characteristic  frequency  doublets  which  were  split  by  about 


half  the  driving  frequency  (f^)  such  that  f-^  +  f2  =  f^  were 
observed.  [20,  21].  Adler  and  Breazeale  [22]  used  Mathieu's 
equation  as  applied  to  the  case  of  a  cavity  with  a  moving 
boundary  to  arrive  at  a  threshold  value  for  subharmonic 
generation.  This  threshold,  expressed  in  terms  of  physical 
parameters,  is 


h  >  ca/u) 

where 

,  _  _A_  _  peak  boundary  amplitude 

Lq  mean  length  of  cavity 

c  =  speed  of  sound  in  the  medium 

a  =  total  absorption  coefficient 

a)  =  0)^/2  =  subharmonic  angular  frequency. 

Yen  [23] 

used  an  acoustic  interferometer  driven  at  1.5  MHz 

with  water  as  the  medium  to  investigate  the  generation  of  a 
subharmonic  doublet  in  two  lower  order  modes.  His  threshold 
for  subharmonic  generation  is 


where 

p  .  8  PC2  ,  1  „  51  S2  .  1/2 

TH  1  +  B/A  Q2  O)01  0)02J 

=  threshold  pressure  for  subharmonic 

generation  in  a  one- dimens ional 

closed  region 

p  =  density  of  the  fluid 

c  =  speed  of  sound 

B/A  =  parameter  of  nonlinearity. 
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The  Q's  are  the  quality  factors  associated  with  the  subharmonic 
modes.  The  6's  are  detuning  factors  such  that 

“01  ■  “l  +  S1 
“02  '  “2  *  S2 
“03  =  “3  *  {3 

where  and  are  the  angular  frequencies  of  the  sub¬ 

harmonic  doublet,  (jOq^  is  the.  driving  frequency  and  the  co^'s 
are  the  angular  frequencies  of  the  resonant  modes. 

It  is  assumed  that  although  the  application  of  Yen's 
threshold  to  the  case  in  which  a  single  subharmonic  component 
is  generated  in  air  at  exactly  half  the  driving  frequency 
may  not  be  entirely  valid,  the  error  introduced  is  not  more 
than  a  factor  of  about  two.  By  also  assuming  that  the  detun¬ 
ing  factors  are  negligibly  small,  =  C^,  B/A  =  2(3-1),  and 
Yen's  threshold  can  be  expressed  as 


l  =  8PC2  r  1  -I 

TH  Q1  L  23-1  J 

■  -^t6-86> 


where  3  =  1.2  in  air.  Recognizing  that  the  Mach  number  of 
the  driven  mode  is 

. .  _  U2  _  P2 

^2  c  'pcT  ’ 

the  threshold  can  be  expressed  as 


M^Q-l  =  6.86 
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The  quantity  characterizes  the  degree  of  the  finite 

amplitude  interaction  and  is  known  as  the  strength  parameter. 
Expression  of  Yen's  threshold  in  this  form  allows  for  direct 
comparison  to  the  theoretically  predicted  threshold  of 
Coppens  [24]  which  is  based  on  similar  mechanisms. 

Coppens '  model  begins  with  the  solution  to  the  one 
dimensional  nonlinear  wave  equation.  It  is  assumed  that  the 
harmonic  components  of  the  distorted  driven  waveform  are  not 

coupled  with  the  subharmonic  and  are  therefore  negligible. 

\ 

The  model  predicts  that  above  a  certain  threshold  the  driving 
level  will  be  sufficient  to  stimulate  a  subharmonic  response 
at  exactly  half  the  driving  frequency.  This  threshold  occurs 
when  the  strength  parameter  is  2.0. 

This  investigation  was  conducted  with  the  objective  to 
experimentally  test  the  theory  of  Coppens  and  others  in  the 
area  of  finite  amplitude  generation  of  subharmonic  frequen¬ 
cies  . 
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III.  EXPERIMENTAL  CONSIDERATIONS 


A.  APPARATUS 

A  block  diagram  of  the  experimental  system  is  shown  in 
Figure  1.  The  resonant  tube,  containing  air  at  ambient 
pressure  and  temperature,  is  six  ft  (1.83  m)  in  length 
with  an  inside  diameter  of  2.250  in.  (5.72  cm)  and  has  steel 
walls  1.125  in.  (2.86  cm)  thick. 

Acoustic  oscillations  were  excited  by  means  of  a  piston 
motion  at  one  end  of  the  tube.  Two  pistons  were  used  during 
the  course  of  the  investigation.  Initially,  a  piston  having 
a  total  mass  of  0.528  kg  and  an  internally  mounted  Endevco 
Model  2215  accelerometer  was  used.  The  output  voltage  of 
the  accelerometer  was  measured  on  a  Hewlett  Packard  Model 
400D  Vacuum  Tube  Voltmeter.  The  output  of  the  voltmeter  was 
displayed  on  a  Tectronic  Model  565  Dual-Beam  Oscilloscope. 
Piston  accelerations  of  up  to  42  g  were  obtained  in  this 
configuration.  To  achieve  the  highest  possible  piston  accel¬ 
erations,  a  light-weight  piston  having  a  mass  of  0.073  kg 
was  subsequently  used.  The  accelerometer  could  not  be 
employed  in  the  light-piston  configuration.  However,  travel¬ 
ing  microscope  observations  of  the  piston  in  open  air  indicated 
that  peak  accelerations  up  to  92  g  were  obtainable.  Each 
piston  had  a  single  Vaseline- lubricated  0-ring  recessed  into 
the  piston  wall  to  enable  a  good  seal  with  the  tube. 
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Figure  1.  Apparatus  Block  Diagram 


The  piston  was  driven  by  an  M-B  Electronics  Model  1500 
Exciter  powered  by  two  M-B  Electronics  Model  2120  MB  power 
amplifiers  connected  in  parallel.  The  driving  signal  was 
generated  with  a  General  Radio  1161-A  Coherent  Decade  Fre¬ 
quency  Synthesizer.  With  this  synthesizer  it  was  possible 
to  select  a  frequency  to  within  0.001  Hertz. 

The  pressure  waveform  within  the  tube  was  sensed  by  two 
condenser  microphones.  At  the  end  of  the  tube  an  Altec  21-BR 
microphone  having  a  diameter  of  one-half  inch  (1.27  cm)  was 
threaded  into  a  1.25  inch  thick  aluminum  end  cap  with  the 
diaphragm  of  the  microphone  flush  with  the  inside  face  of  the 
end  cap.  The  end  cap  was  tightly  bolted  to  the  end  of  the 
tube  with  five  studs.  This  Altec  microphone  was  used 
primarily  to  determine  the  sound  pressure  level  within  the 
tube . 

Because  the  apparatus  was  to  be  driven  at  the  second 
harmonic  of  the  tube  while  looking  for  a  nonlinear  generated 
first  harmonic  expected  to  be  very  much  weaker  than  the  signal 
at  the  driving  frequency,  a  second  microphone  was  mounted  at 
a  point  where  the  standing  pressure  wave  of  the  driven  mode 
has  a  node.  There  are  two  such  points;  one  at  L/4  and  the 
other  at  3L/4  where  L  is  taken  to  be  the  distance  from  the 
piston  face  to  the  end  cap.  (See  Figure  2).  At  these  points 
the  magnitude  of  the  subharmonic  is  0.707  times  its  peak 
value.  The  second  microphone  was  therefore  mounted  at  the 
optimum  point  which  was  furthest  from  the  driver.  A  Bruel 
and  Kjaer  Type  4136  Microphone  was  chosen  for  its  small 
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Figure  2.  Standing  Wave  Diagram 
for  Resonant  Tube 
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diameter  (only  1/4  inch)  and  good  frequency  response  (flat 
to  within  ±0.5  dB  from  80  Hz  to  20  KHz).  The  small  diameter 
was  helpful  in  reducing  the  irregularity  produced  by  mounting 
a  microphone  with  a  flat  diaphragm  in  the  curved  wall  of  the 
tube.  In  order  to  obtain  a  tight  seal  between  the  microphone 
and  the  tube,  an  aluminum  casing  was  constructed  to  securely 
house  the  microphone  with  0-ring  seals.  The  casing  was 
tightly  threaded  into  the  tube  so  that  the  center  of  the 
microphone  diaphragm  was  flush  with  the  inside  wall.  This 
placement  of  the  microphone  reduced  the  signal  at  the  driving 
frequency  by  40  dB. 

After  their  respective  preamplifiers  (see  Figure  1) ,  the 
two  microphone  output  voltages  were  monitored  on  a  pair  of 
Hewlett  Packard  Model  400D  Vacuum  Tube  Voltmeters  and  on  the 
second  and  third  inputs  of  the  Tektronic  Model  565  Dual  Beam 
Oscilloscope.  Each  microphone  was  also  connected  to  a  Hewlett 
Packard  Model  302A  Wave  Analyzer.  Whenever  a  precise  frequency 
analysis  of  the  Bruel  and  Kjaer  wall-mounted  microphone  was 
taken,  a  Ceneco  80593  oscillator  was  connected  to  a  Hewlett 
Packard  Model  521C  Frequency  Counter  and  to  the  wave  analyzer 
which  was  to  be  used  in  the  analysis.  The  oscillator  was 
used  to  align  the  filter  to  the  desired  frequency  before  each 
measurement  and  was  disconnected  while  the  data  was  being 
taken.  (See  section  on  run  procedures.) 

A  Fairchild  Model  766H  Dual-Channel  Oscilloscope  was  used 
to  monitor  the  filtered  output  of  each  wave  analyzer.  A 
Clevite  Model  220  strip  recorder  was  connected  to  the  recorder 
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output  of  each  filter  to  obtain  a  permanent  record  of  the 
data  run. 

B.  CALIBRATION  OF  EQUIPMENT 

Determination  of  the  strength  parameter  and  the  relative 
strength  of  the  subharmonic  required  the  absolute  calibration 
of  both  microphone  systems.  The  calibration  of  the  acceler¬ 
ometer  was  useful  in  determining  the  theoretical  strength 
parameter  achievable  by  the  system. 

Calibration  of  the  microphone  systems  included  the  cali¬ 
bration  of  the  microphones,  their  respective  preamplifiers, 
and  the  VTVM's.  The  Bruel  and  Kjaer  microphone  system  was 
calibrated  using  a  Bruel  and  Kjaer  Model  4240  Pistonphone 
to  produce  a  known  acoustic  pressure  level  of  (124  ±  0.2)  dB 
re  0.0002 y  bar  at  a  frequency  of  250  Hz.  The  sensitivity 
was  found  to  be 

=  (-76.1  ±  0.2)  dB  re  1  volt/y  b. 

=  (1.56  ±  0.04)  x  10'3  volt/Nt  /m2) . 

The  Altec  microphone  was  fitted  with  three  different 
condenser  heads  during  the  course  of  the  investigation. 
Calibration  of  the  Altec  system  was  performed  with  a  modified 
General  Radio  Model  1307-A  Transistor  Oscillator  which  pro¬ 
duces  a  known  acoustic  pressure  level  of  (120  ±  1)  dB  re 
0.0002  y  b  at  a  frequency  of  400  Hz.  The  sensitivities  for 
the  Altec  21  BR  microphone  with  the  three  condenser  heads 
were  found  to  be: 
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21  BR  150:  SM  =  (-62  ±  1) dB  re  lV/yb 

=  (7.9  ±  0.9)  x  10'3  V/(Nt/m2) 

21  BR  200-3:  SM  =  (-87  ±  l)dB  re  lV/yb 

=  (4.5  ±  0.5)  x  10'4  V (Nt/m2 ) 

21  BR  220-3:  SM  =  (-102  ±  l)dB  re  lV/yb 

=  (7.9  ±  0.9)  x  10'5  V/ (Nt/m2 ) 

As  a  check  of  the  relative  calibrations  of  the  two  micro¬ 
phones,  the  tube  was  excited  in  its  third  overtone  at  approx¬ 
imately  378  Hz.  As  shown  in  Figure  2,  the  third  overtone 
has  a  pressure  antinode  at  both  microphones  so  that  the 
acoustic  pressures  measured  by  each  microphone  should  be 
equal.  It  was  found  that  the  two  microphones  gave  acoustic 
pressure  levels  that  were  within  0.2  dB  of  being  equal.  This 
is  well  within  the  calibration  uncertainty. 

The  accelerometer  was  calibrated  by  direct  measurement  of 
the  displacement  amplitude  of  the  piston  face  using  a  travel¬ 
ing  microscope  and  a  strobe.  Ten  measurements  at  a  frequency 

of  190.000  Hz  and  an  accelerometer  output  of  0.120  volts  (rms) 

-4 

resulted  in  an  average  peak  displacement  of  (3.49  ±  0.03)  x  10  m 
and  a  sensitivity  of 

SA  =  (6.82  ±  0.05)  x  10'4  V/ (m/sec2) . 

The  traveling  microscope  and  strobe  were  also  used  to 
measure  the  maximum  obtainable  acceleration  of  the  lighter 
piston.  While  driving  in  open  air  at  a  frequency  of  190.000  Hz 
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a  peak  to  peak  displacement  amplitude  of  0.126  cm  was  measured. 
This  corresponds  to  a  peak  acceleration  of  91.6  g. 

C.  SYSTEM  ALIGNMENT 

Because  the  amount  of  subharmonic  was  expected  to  be  very 
small,  it  was  necessary  to  carefully  align  the  piston  before 
each  run.  The  system  was  designed  so  that  the  resonant  tube 
was  fixed  while  the  driver  assembly  was  mounted  on  a  cradle 
with  elevation  screws  that  allowed  for  three-dimensional 
adjustment.  The  entire  cradle  assembly  rested  on  a  5/16  inch 
pad  of  stiff  acoustic  rubber  to  reduce  creep  at  high  driving 
levels.  The  resonant  tube  was  supported  on  a  separate  table 
by  a  pair  of  foam  rubber  lined  wooden  blocks. 

The  alignment  procedure  consisted  of  exciting  the  tube 
at  a  low  level  while  monitoring  the  accelerometer  waveform 
on  the  dual  beam  oscilloscope.  A  sine  wave  from  the  synthe¬ 
sizer  was  displayed  on  a  separate  trace.  The  cradle  elevation 
screws  were  adjusted  until  the  two  waveforms  appeared  to  be 
identical.  Measurements  of  the  harmonic  distortion  indicated 
that  this  visual  alignment  technique  yields  an  alignment  for 
which  the  second  harmonic  distortion  in  the  piston  motion  is 
less  than  0.3%  of  the  fundamental  and  the  third  harmonic 
distortion  is  less  than  0.5%  of  the  fundamental.  These  values 
are  comparable  to  those  obtained  by  Beech  and  Winn  for  a  sim¬ 
ilar  apparatus.  Any  subharmonic  distortion  was  below  the 
noise  level  of  the  accelerometer  and  was  measured  to  be  less 
than  0.05%  of  the  fundamental. 
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Because  the  accelerometer  was  not  used  in  the  light 
piston  configuration,  alignment  could  only  be  accomplished 
by  adjusting  the  cradle  elevation  screws  until  the  tone 
produced  by  the  piston  sounded  pure.  Measurements  which 
were  taken  while  using  the  heavier  piston  and  accelerometer 
indicated  that  alignment  by  ear  alone  results  in  a  piston 
motion  with  less  than  0.5%  second  harmonic  distortion  and 
less  than  1.0%  third  harmonic  distortion.  Subharmonic  dis¬ 
tortion  is  less  than  0.05%  of  the  fundamental. 

D.  MICROPHONE  OUTPUT  AND  STRENGTH  PARAMETER 
The  strength  parameter  is  defined  as 


SP  =  M2  6 Qx 


where 

M2  =  P2/pc2 


is  the  Mach  number  of  the  driven  mode,  a  quantity  which 
characterizes  the  degree  of  the  finite  amplitude  interaction, 
Q-^  is  the  quality  factor  of  the  subharmonic  mode  and  B  is  1.2 
for  air.  the  acoustic  pressure  of  the  driven  mode,  can 

be  calculated  from  the  microphone  output  if  the  microphone 
sensitivity,  Sn,  is  known.  Using  the  relationships 


S  =  V  /P, 
n  n  2  rms 


P?  =  P,/  ST 
2  rms  2 


the  strength  parameter  may  be  rewritten  as 


SP  =  (1.19  x  10"5)  V  Q./S 
v  J  n  X1  n 


23 


where  Vn  is  the  rms  voltage  of  the  microphone  filtered  at 
the  driving  frequency  and  Sn  is  expressed  in  units  of 
volt/ (NT/m2 ) . 

The  strength  parameter  can  also  be  calculated  from  the 
piston  acceleration  (for  infinitesimal  waves  at  resonance) 
since 


M2  =  AQ/ai'aLc  , 

where 

a  =  Trf/Q2  c 

is  the  absorption  coefficient  for  the  driven  mode,  L  is  the 
tube  length,  and  AQ  is  the  acceleration  amplitude  of  the 
piston.  Aq  is  determined  from  the  accelerometer  output 
voltage  and  sensitivity  from  the  relation 

Ao  -  VSA 

where 

VA  -  ^  VA(rms )• 

The  theoretically  predicted  Mach  number  is  therefore 

M-  =  57.4  VA,  . Q„/f 2 
2  A(rms)x2 

where  V.  ,  ..  is  the  rms  voltage  output  of  the  accelerometer 

A (rms)  0  r 

and  f  is  the  driving  frequency.  This  gives  a  theoretical 
strength  parameter  of 

SP  ■  68‘9  VA(ms)Wf2 

which  should  be  valid  as  long  as  finite  amplitude  effects 
are  negligible. 

All  quality  factors  were  determined  from  the  relation 

0  =  fu  *  f2 

2  (fu  “  f  2) 
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where  f  and  f~  are  the  frequencies  of  the  half  power  points 
11  z 

above  and  below  resonance. 
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IV.  DATA  COLLECTION  PROCEDURES 


A.  PRERUN  PROCEDURES 

Beech  reported  that  the  resonant  frequency  of  the  tube 
changed  by  0.16  Hz  for  each  degree  centigrade  of  temperature 
variation.  Winn  found  that  by  allowing  a  warm-up  period  of 
several  hours  and  by  taking  his  data  in  the  evening  when  room 
temperature  was  relatively  constant,  the  resonant  frequency 
drift  was  of  the  order  of  0.01  Hz  per  hour.  Measurements  for 
this  investigation  were  taken  in  the  evening  following  a  six 
hour  or  longer  warm-up  period  in  which  all  equipment  was 
energized  and  the  tube  was  excited  at  its  second  harmonic. 
Resonant  frequency  drifts  were  found  to  be  consistent  with 
Winn's  data  and  were  considered  insignificant  over  the  dur¬ 
ation  of  the  data  run. 

Piston  alignment  was  accomplished  during  the  warm-up 
period  as  explained  in  Section  III.C.  The  wave  analyzers 
were  calibrated  to  agree  with  the  VTVM's  for  the  same  pure 
tone  input. 

Immediately  preceding  the  run  sequence  the  tube  was 
excited  at  very  low  levels  such  that  the  strength  parameter 
was  less  than  0.20  while  the  half  power  frequencies  needed 
to  determine  the  quality  factors  and  absorption  coefficients 
of  the  first  and  second  harmonics  were  measured. 
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B.  RUN  SEQUENCE 


The  run  sequence  involved  the  collection  of  data  which 
were  necessary  to  obtain  a  spectral  analysis  of  the  Bruel  and 
Kjaer  microphone  output  while  the  tube  was  excited  at  various 
strength  parameters  for  frequencies  near  the  resonance  fre¬ 
quency  of  the  first  overtone  of  the  tube.  The  first  wave 
analyzer,  in  its  Automatic  Frequency  Control  (AFC)  mode,  was 
used  to  measure  the  Altec  microphone  output  at  the  driving 
frequency.  This  allowed  determination  of  the  strength  par¬ 
ameter  of  the  driven  mode.  A  second  wave  analyzer,  in  its 
normal  mode,  was  used  to  measure  the  frequency  spectrum  of 
the  Bruel  and  Kjaer  microphone.  Data  points  were  taken  in 
five  Hertz  increments  from  30  to  220  Hz.  Precise  selection 
of  the  center  frequency  of  the  seven  Hertz  bandpass  filter 
was  accomplished  by  maximizing  the  filter  response  to  a 
calibration  tone  produced  by  the  Ceneco  oscillator  (see 
Figure  1) .  The  filter  was  aligned  in  this  manner  for  each 
frequency  increment.  he  oscillator  output  was  disconnected 
from  the  filter  while  data  were  taken.  At  each  frequency 
increment  the  voltage  level  of  the  Bruel  and  Kjaer  microphone 
was  recorded  from  the  wave  analyzer  for  various  strength 
parameters.  Near  the  frequency  of  the  subharmonic,  where  the 
voltage  level  was  quite  variable,  a  strip  recorder  was  con¬ 
nected  to  the  wave  analyzer  recorder  output  to  produce  a 
permanent  record  of  the  fluctuations  as  a  function  of  time. 
The  strip  recordings  were  later  analyzed  to  determine  average 
and  peak  voltage  levels. 
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V.  RESULTS 


Table  I  lists  characteristic  values  for  resonance  fre¬ 
quency  and  the  associated  quality  factor  and  absorption 
coefficient  for  the  resonant  tube  at  various  modes  in  both 
the  heavy-  and  light-piston  configurations.  Lower  values  of 
a,  which  were  determined  by  Beech  and  Lane  for  the  same  tube, 
were  duplicated  by  replacing  the  1.25  inch  thick  end  cap  with 
the  two-inch  thick  cap  which  was  used  in  earlier  investiga¬ 
tions  to  house  the  Bruel  and  Kj aer  microphone. 

Figure  3  illustrates  the  relationship  between  strength 
parameter  and  accelerometer  output  voltage  as  predicted  from 
first  order  theory  and  as  measured  with  the  microphone.  It 
is  important  to  note  the  increasing  divergence  between  the 
predicted  and  measured  curves  as  the  piston  acceleration  is 
increased.  The  decrease  in  slope  of  the  measured  curve 
represents  the  progressive  increase  of  energy  contained  in 
the  harmonic  components  of  the  distorted  waveform.  These 
components  are  alternately  in  and  out  of  phase  with  the 
driving  frequency,  thereby  resulting  in  a  decrease  in  pressure 
at  the  microphone  from  what  would  be  present  if  all  the  energy 
were  in  the  driven  mode.  Figure  3  also  illustrates  that  in 
the  heavy-piston  configuration  strength  parameters  greater 
than  2.0  cannot  be  obtained.  Greater  piston  accelerations 
were  achieved  using  the  lighter  piston  and  a  maximum  strength 
parameter  of  2.89  was  measured.  This  exceeds  the  threshold 
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for  subharmonic  generation  as  postulated  by  Coppens . 

(SP  _>  2.0).  However,  it  is  considerably  below  the  threshold 
which  was  extracted  from  the  theory  of  Yen.  (SP  >  6.86). 
Since  the  peak  displacement  amplitude  of  the  light  piston  was 
measured  to  be  0.063  cm  at  the  maximum  driving  level,  the 
apparatus  is  operating  below  the  threshold  put  forth  by  Adler 
and  Breazeale  which  requires 


where  =  peak  displacement  amplitude  of  piston 

Lq  =  tube  length 

=  angular  frequency  of  subharmonic. 

For  this  apparatus  Adler  and  Breazeale 's  threshold  reduces 
to  XQ  >  1.92  cm. 

Figures  4  through  10  illustrate  the  frequency  spectrum 
of  the  Bruel  and  Kj aer  microphone  output  for  various  strength 
parameters  ranging  from  zero  to  2.89.  The  three  major  features 
of  Figures  5  through  10  include  the  large  peak  at  the  driving 
frequency,  a  much  smaller  and  variable  peak  at  half  the  driv¬ 
ing  frequency,  and  a  curious  peak  at  70  Hz.  This  70  Hz  peak 
is  thought  to  be  related  to  60  Hz  noise  and  in  Figure  4 
(strength  parameter  =  0)  it  does  occur  at  60  Hz.  The  mechanism 

for  the  translation  of  this  peak  to  70  Hz  for  strength  param¬ 

eters  greater  than  zero  is  not  known. 

The  voltage  level  at  half  the  driving  frequency  varied 
continuously  regardless  of  the  strength  parameter.  All  data 
points  at  95  Hz  indicate  the  average  and  maximum  voltage 
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levels  as  determined  from  strip  recordings  of  the  data.  It 
is  interesting  to  note  from  Figure  4  that  even  when  the  tube 
was  not  being  excited  by  the  piston,  the  peak  at  95  Hz  was 
present  and  quite  variable.  A  possible  explanation  of  this 
observation  is  because  95  Hz  is  the  approximate  frequency  of 
the  first  resonant  mode  of  the  tube,  noise-like  vibrations 
in  the  laboratory  may  have  stimulated  a  weak  response  in  that 
mode . 

Although  the  placement  of  the  Bruel  and  Kjaer  microphone 
resulted  in  a  40  dB  attenuation  of  the  driven  mode,  the  peak 
at  190  Hz  is  by  far  the  dominant  feature  in  Figures  5  through 
10.  Table  II  shows  the  acoustic  pressure  at  95  Hz  as  measured 
by  the  Bruel  and  Kjaer  microphone,  the  acoustic  pressure  at 
the  driving  frequency  as  sensed  by  the  Altec  microphone  and 
the  ratios  of  the  two  pressures  for  each  of  the  strength 
parameters  which  were  investigated.  These  data  do  not  indi¬ 
cate  a  marked  threshold  effect  for  subharmonic  generation 
nor  do  they  indicate  a  linear  relationship  between  the 
magnitude  of  the  response  at  95  Hz  and  the  strength  parameter. 
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Hmw  Pi^tnn  -  ? 

flfl  in  FnrI 

Hole  in 

Tube  Wall  Plugged 

°feeech 

(10Q6n) 

f0(Hz) 

Q 

a (m1 x  100) 

94.904 

73.28 

1.18 

1.1 

189.194 

118.95 

1.46 

1.4 

285.548 

142.85 

1.83 

1.8 

380.977 

163.51 

2.13 

2.1 

476.642 

186.33 

2.34 

2.3 

572.338 

203.03 

2.57 

2.8 

Heavy  Piston  -  1.25”  End  Cap 

Q 

a  (m-1  x  100) 

94.720 

73.48 

1.18 

189.831 

114.22 

1.52 

284.964 

•137.93 

1.89 

380.209 

158.03 

2.20 

475.560 

181.72 

2.39 

571.004 

200.63 

2.60 

Light  Piston  -  1.25  in. 

End  Cap 

f0(Hz) 

Q 

a  (rr 1x  100) 

94.410 

56.10 

1.54 

189.252 

96.05 

1.80 

Table  I.  Resonant  Frequencies,  Quality  Factors,  and 
Absorption  Coefficients  for  Resonant  Tube 


51 


JsasuiBJBj  q^Susjas 


32 


0.00  0.02  0.04  0.06  0.08  0.10  0.12  0.14  0.16  0.18  0.20 

Accelerometer  Voltage  (Vnns) 

Pigure  3.  Strength  Parameter  vs.  Accelerometer  Voltage 
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Figure  4.  Frequency  Spectrum  of  Bruel  and  Kjaer 
microphone . 


Strength  parameter  -  0.707 

Heavy  Piston  Configuration 
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;igure  5.  Frequency  Snectrum  of  Bruel  and  Kjaer  Microphone. 
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Figure  6.  Frequency  Spectrum  of  Bruel  and  Kjaer  Microphone. 
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Figure  7.  Frequency  Spectrum  of  Bruel  and  Kjaer  Microphone. 
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Figure  8.  Frequency  Spectrum  of  Bruel  and  Kjaer  Microphone 


33 


Figure  9.  Frequency  Spectrum^ f  Bruel  and  Kjaer  Microphone 
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Figure  10.  Frequency  Spectrum  of  Bruel  and  Kjaer  Microphone 


Table  II.  Comparison  of  Acoustic  Pressures  at  190  Hz 
and  95  Hz  for  Various  Strength  Parameters 


Strength 

Parameter 

P190  (Nt//m2) 

Pg5  (Nt/m2) 

(dB)  p95 

20  log  (p^-) 

*190 

HEAVY  PISr 

rON  CONFIGURATION 

0.707 

1.14  x  103 

4.43  x  10'2 

-88.2 

1.18 

1.91  x  103 

0.125 

-83.7 

1.55 

2.51  x  103 

9.05  x  10‘2 

-88.8 

LIGHT  PIS! 

rON  CONFIGURATION 

1.04 

2.21  x  103 

0.126 

-84.9 

2.09 

4.43  x  103 

0.157 

-89.0 

2.89 

6.13  x  103 

0.433 

-83.0 

P190  =  Pressure  at  190  Hz  as  computed  from  filtered  output  of  Altec 
microphone . 

P95  =  AveraSe  pressure  at  95  Hz  as  computed  from  filtered  output  of 

Bruer  and  Kjaer  microphone. 
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VI.  CONCLUSIONS 


Detailed  measurements  of  the  pressure  waveform  within 
the  resonant  tube  indicate  that  strength  parameters  up  to 
2.89  were  obtained.  Although  the  driving  levels  were  not 
high  enough  to  test  the  theories  of  Adler  and  Breazeale  and 
Yen,  the  theoretical  threshold  for  subharmonic  generation, 
as  developed  by  Coppens,  was  exceeded.  Experimental  obser¬ 
vations  indicate  that  there  is  no  noticeable  threshold  effect 
in  the  range  of  strength  parameters  achievable  by  the  system. 
Therefore,  Coppens'  theory  was  not  verified.  Any  subharmonic 
which  may  have  been  present  was  more  than  83  dB  below  the 
acoustic  pressure  level  of  the  driven  mode. 
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